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We compare the roles of fidelity and entanglement in characterizing renormalization group fiows 
and quantum phase transitions. It turns out that the scaling parameter extracted from fidelity 
for different ground states succeeds to capture nontrivial information including stable and unstable 
fixed points, whereas the von Neumann entropy as a bipartite entanglement measure (or equivalently, 
majorization relations satisfied by the spectra of the reduced density matrix along renormalization 
group fiows) often fails, as far as the intrinsic irreversibility-information loss along renormalization 
group flows-is concerned. We also clarify an intimate connection between the von Neumman entropy, 
majorization relations, and fidelity. The relevance to Zamolodchikov's c theorem is indicated. 
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Quantum phase transitions (QPTs) [H are a topic of 
current interest, subject to intense study in condensed 
matter systems. Conventionally, a QPT is described in 
terms of a Hamiltonian and its spectrum. Most examples 
of the well-studied QPTs fit into the Landau-Ginzburg- 
Wilson paradigm, with the central concept being a local 
order parameter. The non-zero expectation value of a lo- 
cal order parameter characterizes a symmetry-breaking 
phase that only occurs for a system with infinite num- 
ber of degrees of freedom, in contrast to QPTs result- 
ing from a level crossing in a finite size system. How- 
ever, there exist phases beyond symmetry-breaking or- 
ders, which yields exotic QPTs characterized by topo- 
logical/quantum orders [1]. Unconventional QPTs also 
occur in matrix product systems y. 

Another picture emerges due to latest advances in 
quantum information science, which attempts to char- 
acterize QPTs from both entanglement \^ M , S 0, Q 
(for a review, see (al) and fidelity ^ 

perspec- 
tives. Remarkably, for quantum spin chains, the von 
Neumann entropy, as a bipartite entanglement measure, 
exhibits qualitatively different behaviors at and off crit- 
icality Q. On the other hand, it has been shown that 
fidelity may be used to characterize QPTs, regardless of 
what type of internal order is present in quantum many- 
body states Uj . A basic point is that there always exists 
an order parameter (local/nonlocal) in a system under- 
going a QPT. In principle, such an order parameter may 
be constructed systematically [l^]. This results in the 
orthogonality of different ground states, due to state dis- 
tinguishability [l^. Thus, it is justified that the scaling 
parameter may be extracted from fidelity to character- 
ize QPTs. In fact, nontrivial information about QPTs 
including stable and unstable fixed points along renor- 
malization group (RG) flows may be revealed solely from 
ground states, and critical exponents can be extracted 
by performing a finite size scaling analysis from such a 
scaling parameter [3]. 



An intriguing question is to clarify the connection be- 
tween fidelity and entanglement approaches to QPTs. An 
apparent difference between entanglement measures and 
fidelity is that entanglement measures involve partitions, 
whereas the fidelity approach treats systems as a whole. 
Therefore, no matter what bipartite entanglement mea- 
sures are used, some physical information will be lost, 
since the whole is not simply the sum of the parts. On 
the other hand, one may expect that the scaling param- 
eter extracted from fidelity should capture information 
that is lost in bipartite entanglement measures. 

In this paper, we compare the roles of fidelity and en- 
tanglement in characterizing RG flows and QPTs. It is 
found that the scaling parameter extracted from fidelity 
for different ground states succeeds to capture nontrivial 
information including stable and unstable fixed points, 
whereas the von Neumann entropy as a bipartite entan- 
glement measure (or equivalently, majorization relations 
satisfied by the spectra of the reduced density matrix 
along RG fiows) often fails, as far as the intrinsic irre- 
versibility, i.e., information loss along RG flows is con- 
cerned. We also clarify an intimate connection between 
the von Neumman entropy, majorization relations and 
fidelity. The relevance to Zamolodchikov's c theorem is 
indicated, if a system is conformally invariant at transi- 
tion points,. 

Generalities. For a quantum spin chain described by a 
Hamiltonian H{X), with A a control parameter, it is well- 
established J^] that its ground state may be rep- 
resented in terms of the so-called matrix product states 
(MPS) [Tt^. Suppose the system is translationally invari- 
ant, then takes the form, 



|V(A)) =Tr(A,,(A).--A,JA))|*i-..Zi), 



(1) 



where {Ai} is a set of dD x D matrices, with d being the 
dimension of the local Hilbert space at each lattice site, 
and D the dimension of the bonds in the valence bond 
picture The fidelity -F(A,A') = \{ip{X')\ip{X))\ for 
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two ground states |^/'(A)) and IV'(-A')) corresponding to 
different values of the control parameter A is 



i^(A,A') ^ Tr(i;^(A,A')), 



(2) 



with the (generalized) transfer matrix £'(A, A') = 
Es=i^''*(^') ® ^'{^)- Therefore, it is straightforward 
to evaluate the fidelity F by diagonalizing the transfer 
matrix i?(A, A'). For large enough L's, the fidelity F 
scales as F(A, A') ^ d^(A, A'), where d(A, A') is the scal- 
ing parameter. In the thermodynamic limit, (i(A, A') may 
be defined as 



lnd(A, A') = lim lni^(A, A')/i. 



(3) 



In fact, the scaling parameter (i(A, A) is the largest eigen- 
value of the transfer matrix E{X,X'), if there is no pair 
of eigenvalues that are complex conjugate each other 
(otherwise, a fast oscillating part should first be fac- 
tored out [11]). Obviously, it satisfies the properties in- 
herited from the fidelity -F(A,A'): (a) d{X, X) = 1; (b) 
d(A, A') = d(A', A); and (c) < d(A, A') < 1. 

Although the system is in a pure state, a subsystem A 
consisting of N consecutive spins is generically in a mixed 
state described by the reduced density matrix pa{X) — 
Ti bp{X), with the density matrix p{X) — |'(/'(A))('0(A)|, 
and the trace is taken over the subsystem B comple- 
mentary to the subsystem A. The fidelity Fa (A, A') for 
two reduced density matrices Pa{X) and Pa{X') is de- 
fined as Fa (A, A') = Ti^J pK^{X)pa{X')p](^{X). For large 
enough L's and iV's, with a fixed N/L [l^. Fa scales 
as Fa (A, A') ~ d^{X,X'), where dA(A, A') is the scaling 
parameter for the subsystem A. Formally, in the ther- 
modynamic limit, dA{X,X') may be defined as 



lndA(A,A')= lim In Fa (A, A')/iV. 



N 



(4) 



The scaling parameter dA{X, X) also enjoys the proper- 
ties: (a) dA(A,A) = 1; (b) dA(A, A') = dA(A',A); and (c) 
< dA{X,X') < 1. The scaling parameters d(X,X') and 
dA(A, A') satisfy 



dA{X,X')>d{X,X'). 



(5) 



This follows from the well-known inequality: Fa(A, A') > 
F(A, A') [l^. Although it is straightforward to cal- 
culate d{X, A') in the MPS representation, it is a 
formidable task to evaluate dA{X,X'). This is due to 
the fact that dA(A, A') depends not only on the spec- 
tra but also on the eigenvectors of the reduced den- 
sity matrices pa(A) and pa(A'). From the Schmidt de- 
composition, one may expect that dA{X,X') takes the 
form dA(A, A') = Y.a,/3 9a(} (A, X')yjuja (A) y^w^j ( A') , with 
a,/9 = 1, • • • , F* and {wa{X)} the (normalized and non- 
increasing ordered) spectra of the reduced density ma- 
trix Pa, and gaf3{X,X') a model-dependent real ten- 
sor, satisfying gct^(A, A') = gp^ai^' , X). The difficulty 



arises from the dependence of ga/3{X,X') on A and A'. 
Thus, instead of dA(A, A'), we may seek a quantity 
dE{X,X') = Y.a,f3 ^ai3 (A) ^cj^ ( A') , which is uni- 
versal in the sense that the real tensor F^^ is inde- 
pendent of A and A'. The only choice is Eajs — Sap, 
if one requires that d£;(A, A') enjoys the same proper- 
ties as d(A, A') and dA{X,X'): (a) d£;(A, A) — 1; (b) 
dB(A,A') dE{X',X)- and ( c) < d£;(A,A') < 1. 
Thus dE{X,X') = Y.a V^^^4A)a/I^4A0. We stress that 
dE{X,X') only involves the eigenvalues of the reduced 
density matrix, whereas d{X,X') unveils information en- 
coded in the entire ground state. 

The connection between d{X,X'), dE{X,X') and ma- 
jorization. Suppose a system fiows to two different stable 
fixed points A_ and A-|_ under RG transformations [l^ . 
with an unstable fixed point Ac as a transition point. 
An intriguing fact is that dE^X, A') shares similar behav- 
iors to d(A, A'), if A and A' are in the same phase. To 
demonstrate this, we notice that, generically, there is one 
and only one extreme point for ^^(A, A') if one regards 
dE{X,X') as a function of A for a given A'. That is, when 
A = A', it reaches the maximum 1, as follows from the 
Lagrange multipliers for c?£;(A, A') under the constraints: 
E„Wa(A) = 1 and Ea^a(A') = 1. Therefore, d£(A,A') 
increases with A, until it reaches 1 when A = A', then it 
decreases with A along an RG flow for a fixed A'. 

We stress that the monotonic behaviors of dE{X,X') 
are consistent with majorization relations for the spectra 
of the reduced density matrix along RG flows. It has 
been established that there is a more "fine-grained" 
characterization of entanglement loss along RG flows 
in terms of majorization [2C 

"ife = J2a=l'^a{k = 1,2,- 



2ll . l22l |. which states that 



is non-decreasing along 
RG flows. As shown in Ref. [2^, this simply follows from 
the fact that there is one and only one crossing point a* 
for eigenvalue distributions when one regards LUa as a 
function of the index a [1^. Equivalently, Vac^q > for 
a < a* , and V \uja < for a > a*, where V\ denotes 
the derivative with respect to A along RG flows. Then, 
suppose A and A' are in the same phase and A' is fixed, 

VArfg(A, AO ^ l/2(E a<a> V^„(A')/Wa(A)VAa;a(A) + 

E„>„. V^a(A')/u;a(A)VAa.4A)) > l/2\/ xiE^^^W) ^ 
0, if A flows from Ac to A', because uJa{X')/uJa{X) > 1 
for a < a*, and uia{^')/^^a{X) < 1 for a > a*; 
VArig(A, AO = l/2(E c<a> ^LUaiX')/u;a{X)^xtOaW + 
'Eaya- v/^^a(A')/wa(A)VAu;a(A)) < 0, if A flows from A' 
to a stable fixed point, because uja{X')/u!a{X) < 1 for 
a < a*, and LUa{X')/u!a{X) > 1 for a > a*. 

If A and A' are in different phases, then the mono- 
tonic behaviors of dE{X,X') with A along an RG fiow 
for a fixed A' depends on what type of orders present 
in systems. For our purpose, we discuss the two ex- 
treme cases: = 2 for MPS systems, and D = oo 
for spontaneous symmetry-breaking orders. Since the 
monotonicity is universal in the sense that it should re- 
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main the same for A"s from the same phase, therefore we 
just need to focus on the stable fixed points A' = A±. 
For D — 2, there are only two different choices for 
the reduced density matrix spectra at the stable fixed 
points: {1,0} and {1/2,1/2}, corresponding to unen- 
tangled states and maximally entangled states, respec- 
tively. For the former, V\dE{\X') = \/x^yuJl{X) > 0, 
since a;i(A) is non-decreasing along an RG flow. For 
the latter, VAci^CA, A') = l/^^x{^/^M>^ + V^M>^) = 
l/(2v^)(l/y^J^- 1/V^J;(A))VaWi < 0. For = oo 
as occurs for spontaneous symmetry-breaking orders, 
a remarkable feature is that the reduced density ma- 
trix spectra at critical points tend to vanish, leading to 
the divergence of the von Neumann entropy for a half- 
infinite chain. Thus, De{X,X') tends to when A or 
A' tends to Ac. On the other hand, degeneracies oc- 
cur in the largest eigenvalue of the reduced density ma- 
trix in symmetry-broken phases, and no degeneracy in 
symmetric phase. Therefore, the reduced density ma- 
trix spectra are {1/n, ••• ,1/n, O,---} (n-degeneracies) 
for stable fixed points in symmetry-broken phases and 
{1,0, •••} for stable fixed points in symmetric phases. 
Then, ^^(A, A') as a function of A for a fixed A' mono- 
tonically increases, if A' is in symmetric phases, whereas 
d£;(A, A') increases, until it reaches a maximum, then de- 
creases, if A' is in symmetry-broken phases. 

The fact that c?£;(A, A') shares similar monotonic be- 
haviors to d{X, A') establishes a connection between en- 
tanglement and fidelity approaches to QPTs, since infor- 
mation unveiled in cIe results from the reduced density 
matrix spectra that in turn determine the von Neumann 
entropy E{X) = — J2a '^aW liiuja{X). The latter is non- 
increasing along RG fiows. If a system is conformally 
invariant at transition points, then one may formulate 
an entropic version of Zamolodchikov's c theorem [j^ . 
Therefore, the monotonicity of d{X,X') and dE{X,X') is 
reminiscent of Zamolodchikov's c theorem. 

Quantum XY spin 1/2 chain. The quantum XY spin 
chain is described by the Hamiltonian 



M 



j=-M 



Here crj,^^, and fT| are the Pauli matrices at the j-th 
lattice site. The parameter 7 denotes an anisotropy in 
the nearest-neighbor spin-spin interaction, whereas A is 
an external magnetic field. The Hamiltonian ([6]) may be 
diagonalized exactly (25j . In the thermodynamic limit, 
d{X, A'; 7) takes the form: 



d(A,A';7) = exp[— dalnT{X,X';r,a)], (7) 

where J-{X, A'; 7; a) — cos[d{X; 7; a) — 't?(A'; 7; a)]/2, with 

cos ??(A; 7; a) = (cos a — A)/ ^ (cos a — A)^ + 7^ sin^ a. It 
was shown [ll| that d{X,X';j) exhibits a pinch point at 



(1, 1), i.e., an intersection of two singular lines A = 1 and 
A' = 1. For quantum Ising model in a transverse field 
(7 — 1), all states for A > 1 flow to the product state 
with all spins aligning in the z direction (A = 00), and 
all states for A < 1 flow to the cat state (| ^ • • • ^ 
) + I ^ • • • ^))/V2 (A = 0). We stress that d(A, A'; 7) 
detects the duality between two phases: A > 1 and A < 1, 
because d(A, A'; 1) = d(l/A, 1/A'; 1) [H]. 

The exact spectra of the reduced density matrix for a 
half-infinite chain take the form 12611: 



j(A;7) 



ZiX-n) 



(8) 



where — (2j + l)e if A > 1, and ej = 2je if 
A < 1, and Z(A;7) = (16g/(fc2fc'2))i/24 if A > 1, and 
Z(A;7) = 2(16A:V((7fc'))^/^^ if A < 1, with q = exp(-e). 
Here e = TTK{k')/ K{k), with K{k) being the complete 
elliptic integral of the first kind, k the modulus and 
k' = y/l — fc2. The relation between the modulus k and 
the anisotropy 7 and the transverse field strength A is 
k = 7/VA2 + 72 - 1 if A > 1 and fc = y^J^^T^f^/j ff 
A < 1. Therefore, if A and A' are in the same phase, we 
may derive a closed expression for d_B(A, A'; 7): 



rfij(A,A';7) = 



^(A,V;7) 
V^(A;7)^(A';7)' 



(9) 



where Z(A, A';7) takes the same form as Z(A;7), i.e., 
Z(A,A';7) = (16g/(Pfc'2))i/24 jf a. A' > 1, and 
Z A,A^;7) = 2{16k^/{qk')y^^^ if A, A' < 1, with q = 
^ q{X)q{X'), and k determined from q = exp(— e), with 
e — TTK{k')/K(k). However, if A and A' are in differ- 
ent phases, (i£;(A, A';7) is only available in terms of an 
infinite sum. 

The monotonic behaviors of dsiX, A'; 7) as an example 
for D = 00 with Z2 broken symmetry are numerically 
confirmed. If A and A' are in the same phase, o?£;(A, A'; 7) 
takes values very close to c?(A, A';7) when A and A' are 
away from the critical point. However, in contrast to 
d(A, A'; 7), dsiX, A'; 7) fails to detect the duality between 
two phases for quantum Ising model in a transverse field 
(7 = 1). This is due to the fact that the dual unitary 
transformation connecting the two phases involves non- 
local operations, which do not affect the fidelity (and so 
d(A, A'; 7)), but do change the spectra of the reduced den- 
sity matrix (and so dsiX, X';j)). We note that the von 
Neumann entropy E{X) is decreasing along the RG flows, 
so it succeeds to detect the stable fixed points: A = and 
A = 00, but it fails to detect the duality. 

Quantum spin 1/2 chain with three-body interactions. 
The system is described by a Hamiltonian with three- 
body interactions [27t 



H 



2{g^-l)ataU,~{l+gf<j^,H9~^f<yt'y^.+i<y 



(10) 
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The peculiarity of the model is that it is not conformally 
invariant at the transition point gc = 0: in the thermo- 
dynamic limit, the correlation length diverges, and en- 
ergy gap vanishes, but the ground state energy remains 
smooth 0. As emphasized in Ref. the parameter 
space should be compactified by identifying g — +oo with 
g ~ — oo, due to the fact that H{+oo) = H{—oo). Since 
ground states are an MPS with Ai = (/ — (7^)/2 + (7^ and 
A2 = (/ + (T^)/2 + g(T+ 3], one may extract the scaling 
parameter d as d{g, g') = ^1 + l<?5'l/\/(l + + l^'l) 

if g and g' are in different phases, and d{g,g') — (l + 
VI55'I1Za/(1 + l.9l)(l + WD if 9 and g' are in the same 
phase [nl . From this we read off that there are two tran- 
sition points: g = and 00. All states for positive g flow 
to the product state (5 = 1) with all spins aligning in 
the X direction, and all states for negative g flow to the 
cluster state ^ [g = -1). 

The reduced density matrix spectra for a half-infinite 
chain are identical to those of p, which may be calculate 
exactly by exploiting the freedom in the set {Ai} in order 
to fix the gauge: A,A\ = / and Ap^^ = P 3] . For 
5 > 0, the eigenvalues of the reduced density matrix are 
1/2 ± Y^/(l -I- g), whereas for <? < 0, both of them are 
1/2. Obviously, for positive g, the majorization relations 
are satisfied, since the larger eigenvalue 1/2-1- y^/ (1 -I- g) 
is monotonically increasing when g varies from to 1 , or 
varies from 00 to 1, consistent with the fact that both 
g — Q and g = 00 are unstable fixed points, and g — 
1 is a stable fixed point. However, the trivial spectra 
for negative g fail to identify another stable fixed point 
g = —1, in contrast to d{g,g'). One may check that 
dE{g,g') = d{g,g') if > 0, and dE{g,g') = 1, i.e., 
dsig-g') serves as a trivial upper bound, if g,g' < 0. On 
the other hand, dE{g,g') = l/Vl + .9 if < g < 1 and 
g' < 0, and dE{g,g') = ^/g/{l + g) if 1 < .g < 00 and 
g' < 0. In this case, dE{g,g') is decreasing with g along 
an RG flow for a fixed g' , consistent with the majorization 
relations. The von Neumann entropy E{g < 0) = In 2 
also fails to identify the stable fixed point 5 = — 1. 

Quantum spin 1 chain with two-body interactions. 
The model is a deformation of the celebrated AKLT 
model [l3|, with the Hamiltonian 

H = ;^(2 + g^)SA+, + 2(SiS,+i)2 + 2(4 - g^){S')^ 

■i 

+ [g + 2f{StSU^f + g{g + 2){5,^5f+i, S,S,+i}+. 

(11) 

Note that the AKLT model corresponds to t; = —2. 
The model is very similar to the previous one, although 
the order parameter is non-local. Indeed, the parame- 
ter space should be compactified by identifying g — +00 
with g = —00, since H{-\'Oo) = H{—oo). The ground 
states are an MPS with Ai — — cr^,^2 — and A3 — 
ga+ i, thus dig,g') - ^1 + l55'l/V(l + l5l)(l + l.9'l) 
if g and g' are in different phases, and d{g,g') — (l + 



Vl55'l)/\/(l + l5l)(l + \g'\) if 9 and g' are in the same 
phase. There are two transition points: g = and 00. 
All states for positive g flow to the state with g = 1 , and 
all states for negative g flow to the state with g = —1. 

The reduced density matrix spectra are trivial, since 
p — 1/2 I for both positive and negative g. Therefore, no 
information is unveiled from the spectra (and so the von 
Neumann entropy E{g)), as far as the stable fixed points 
are concerned. Actually, dE{g,g') — 1, and E{g) — In 2. 
The fact that both d and remain the same under 
interchange g ^ —g reflects that the ground states for 
±g are equivalent up to local unitaries 
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